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We study the interaction of a nonlinear spin- wave and magnetic soliton in a uniax- 
ial anisotropic ferromagnet. By means of a reasonable assumption and a straightfor- 
ward Darboux transformation one- and two-soliton solutions in a nonlinear spin-wave 
background are obtained analytically, and their properties are discussed in detail. 
In the background of a nonlinear spin wave the amplitude of the envelope soliton 
has the spatial and temporal period, and soliton can be trapped only in space. The 
amplitude and wave number of spin wave have the different contribution to the 
width, velocity, and the amplitude of soliton solutions, respectively. The envelope 
of solution hold the shape of soliton, and the amplitude of each envelope soliton 
keeps invariability before and after collision which shows the elastic collision of two 
envelope soliton in the background of a nonlinear spin wave. 

PACS numbers: 75.75.+a, 05.45.Yv, 76.50.+g, 73.21.Hb. 

I. INTRODUCTION 

Solitons are the natural results of the nonlinear equation in which the dispersion is com- 
pensated by nonlinear effects. With this compensation solitons can travel over long distances 
with neither attenuation nor change of shape which has the prominent application in high- 
rate telecommunications with optical fibers in the future. For this reason the study of solitons 
has received considerable^ attentions in many fields, such as particle physics, molecular bi- 
ology, nonlinear optics [1'] and condensate physics One of the topics in condensate 
physics is the nonlinear dynamics of local magnetization in Heisenberg spin chain model 
which successfully explains the existence of ferromagnetism and antiferromagnetism at tem- 
peratures below the Curie temperature. Taking into account the spin-spin interactions, the 
nonlinear excitations, such as spin wave and magnetic solitons ^, , are general phenomena 



in ordered ferromagnetic materials. By means of the neutron inelastic scattering [3] and 
electron spin resonance js], the magnetic soliton had already been probed experimentally 
in quasi-one-dimensional magnetic systems. The magnetic soliton, which describes local- 
ized magnetization, is an important excitation in the classical Heisenberg spin chain. In 
particular, the continuum limit for the nonlinear dynamics of magnetization in the classical 
ferromagnet is governed by the Landau-Lifshitz (L-L) equation [9]. This equation governs 
a classical nonlinear dynamically system with novel properties. In a one-dimensional case, 
some types of L-L equation is complete integrable. The isotropic case has been studied in 



various aspects [lOl - 
is also discussed 



mg 



13, 



transformation 



2 L and the construction of soliton solutions of anisotropic L-L equation 
-II 811. In these studies a variety of techniques, such as inverse scatter- 
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Riemann-Hilbert approach 16|, and Darboux transformation 



18| , have been applied to construct the soliton solution of L-L equation. Recently, this 



Heisenberg model has new application in some fields of condensate physics. In rnagnetic 
multilayers the magnetization switching and reversal 19|], domain- wall dynamics 20| and 



magnetic solitons 



21 



22| were discussed analytically in terms of L-L equation with spin- 



torque. Taking into account the light-induced and the magnetic dipole-dipole interactions 



';he ferromagnetic state 



23| are obtained and the nonlinear excitations, such as spin wave 



24j | and magnetic solitons 25|, has been reported in spinor Bose-Einstein condensates in 



deep optical lattice. 

In this paper, we will study the interaction of a nonlinear spin wave and magnetic solitons 
in a uniaxial anisotropic ferromagnet. In terms of a reasonable assumption we transform the 
Landau-Lifshitz equation into an equation of the nonlinear type. By means of a straight- 
forward Darboux transformation we report one- and two-soliton solutions in a nonlinear 
spin-wave background analytically and discuss their properties in detail. 

II. MODEL 



In the classical limit, the dynamics of magnetization of a ferromagnet as a function of 
space and time M [x, t) is determined by the Landau-Lifshitz equation 

dM 2fi. 



dt 



h 



(1) 
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where /Iq is Bohr magneton, the effective magnetic field He// is equal to the variational 
derivative of the magnetic crystal energy with respect to the vector M, 

where the magnetic crystal energy function including the exchange energy, anisotropic energy 
and the Zeeman energy can be written as 

2 

dx 



2 



£x, (3) 



where A is the exchange constant and /3 is the uniaxial anisotropic constant, /3 > corre- 
sponds to easy-axis anisotropy while (3 < corresponds to easy-plane type. Substituting 
Eqs. ^ and Q into Eq. (P, we get 

l~^ = -^Mx ^ - /3M X e3(M ■ e,) - /^qM x B, (4) 

where 63 is the unit vector along the 2;-axis, and B = (0, 0, B). When the magnetic field is 
high enough, the deviation of magnetization from the the direction of the field is small. With 
the consideration that at temperatures well below the Curie temperature the magnitude of 
magnetization is a constant, i.e., the integral of motion = Mq = constant, here Mq is the 
saturation magnetization, we can introduce a single function ^ , instead of two independent 
components of M, 



\1/ = m^; + imy, = y 1 - |vl'|^, (5) 
where m = {mx,niy,mz) = M/Mq. Thus Eq. becomes 

AMo m,— mJ = -(/3Mom, + S)^. (6) 



2/iQ dt \ dx"^ dx"^ 

In the ground state, vector M directs along the anisotropy axis 63. Now we consider 
the small deviations of magnetization from the equilibrium direction which corresponding to 
nv^+niy <ti m^, i.e., <^ 1, then ^ 1 — 1/2 {"^f. In the long-wavelength approximation 
and the case /3 > 0, Eq. (jH]) may be simplified by keeping only the nonlinear terms of the 
order of the magnitude of |?/^|^ ^. As a result, we have the following dimensionless Schrodinger 
equation: 



(7) 
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For convenience we have rescaled space x and time t by 21q and I/cuq, where Iq = ^/ A/jS 
is the characteristic magnetic length and Uq = Pfi^Mo/h is the homogeneous ferromagnetic 
resonance frequency. The spin wave interaction in a ferromagnet with easy-axis anisotropy 
is attractive in nature. This attraction leads to the macroscopic phenomena which are 
associated with the appearance of spatially localized magnetic excitations, i.e., magnetic 
solitons which is admitted in nonlinear equation ([7]). So it is very interesting to investigate 
the dynamics of magnetic soliton in a nonlinear spin wave background which hasn't been well 
explored. To this purpose we have to get soliton solutions of ([7]) embedded in a nonhnear 



26 



has been 



spin-wave background. It should be noted that Darboux transformation 
developed to construct such solutions. The main idea of this method is that it firstly 
transforms the nonlinear equation into the Lax representation, and then in terms of a series 
of reasonable transformations the soliton solutions can be constructed algebraically with an 
obvious seed solution of the nonlinear equation. By employing Ablowitz-Kaup-Newell-Segur 
technique we can construct Lax representation for Eq. ([7]) as follows 



— ^ 
dx 



d_ 

di 



(8) 



where = yiJi 4'2 ) ' superscript "T" denotes the matrix transpose. The lax pairs U 
and V are given in the forms 



U 



AJ + P, 



V = (-iX^ + 02) J- iXP + lt{P^ + -^P)J, 
^ ' 2 ax 



(9) 



with 





^10 ^ 
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B , 
/3Mo' 



where the overbar denotes the complex conjugate. With the natural condition of Eq. ([8]) 
-^ol'ip = ^S^"*/"' i-^-' ~ + ^] — 0' -^l- © can be recovered. Based on the 
Eq. ([8]), we can obtain the general one- and two-soliton solution by using a straightforward 



Darboux transformation 



26|-l28| 
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III. DARBOUX TRANSFORMATION 

In order to clear the correction of soliton solutions we briefly introduce the procedure 
of getting soliton solutions for the developed Darboux transformation. We consider the 
following transformation 

^[1] = (A/-^)^, (10) 
where S = KAK~^, A =diag(Ai, A2) , and i^' is a nonsingular matrix which satisfies 

= JKA + PK. (11) 

Letting ip [1] satisfy the Lax equation 

[1] = f/i^ [1] , (12) 

/ ^1 \ 

where Ui = \ J + Pi, Pi = \ _ , and with the help of Eqs. ([TU]) and ([TT]). we 

obtain the Darboux transformation for Eq. ([7]) from Eq. (fT2!) in the form 

^1 = ^ + 25i2. (13) 

The above equation implies that a new solution of Eq. ([7]) can be obtained if S is known. To 
obtain the expression of 5* it is is easy to verify that, li ip = {^ip^ is a eigenfunction of 

Eq. ([H]) corresponding to the eigenvalue A = Ai, then^ —''P2 i^i ) eigenfunction, 
while with the eigenvalue — Ai. Therefore, K and A can be taken the form 

A-=| |.A=| ■ _ I. (14) 



^2 V^l 

which ensures that Eq. fITT]) is held, we can obtain 




Ssi = -XiSsi + (Ai + Ai) s, / = 1, 2, (15) 



where ijj'^ilj = \'^i\^ + \'^2?' ■ Then Eq. (IT^ becomes 

vl/i = vl/ + 2(Ai + A0^, (16) 

where = (y'lpi 4^2^ to be determined is the eigenfunction of Eq. (jS]) corresponding to the 
eigenvalue Ai for the solution Thus if Eq. ^ is solved we can generate a new solution 
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therefore the new solutions of m from Eq. (j^]), of the Eq. ([7]) from a known solution \l/ 
which is usually called "seed" solution. 

To obtain exact A^-order solution of Eq. , we firstly rewrite the Darboux transformation 
in Eq. (fT6i) as in the form 

where ip [1, A] = {ipi [1, A] , ■ip2 [I5 ^])'^ denotes the eigenfunction of Eq. Q corresponding to 
eigenvalue A. Then repeating above the procedure N times, we can obtain the exact A^-order 
solution 

^. = ^ + 2 y:(A„ + j^^hhMMhM^ (is) 

where 

^[n,X] = {X-S[n-l])---iX-S[l])4j[l,X], 

^ 1,]- \^ , ^y^iM^dZiM 

^si [J\ - -Xjdsi + [Xj + Xj) ^ , 

here ■?/' [j, A] is the eigenfunction corresponding to Xj for ^E'j-i with \l/o = \1/ and s, / = 1, 2, 
j = 1, 2, ■ ■ ■ , n — 1, = 2, 3, ■ ■ ■ , A^. Thus if choosing a "seed" as the basic initial solution, 
by solving linear characteristic equation system ([8]), one can construct a set of new solutions 
of m from Eqs. ([5]), ((ZD and ([l8]). 

In order to get soliton solutions in a nonlinear spin wave background we take the initial 
"seed" solution as 

which corresponding to a nonlinear spin wave solution of m as 

TTLx = Ac cos (kcX—Uct) , 

niy = —Ac sin (kcX—Uct) , 



m, = - (20) 

which satisfies the nonlinear dispersion relation Uc = kc/2 — Ac + 2B/ (/3Mo) + 2, where 
the amplitude Ac is small real constants and kc is the wave number of nonlinear spin wave. 
Substituting ( 1T9|1 into ([7]) and solving the linear equation system ([8]), after the tedious 
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calculation we have the eigenfunction of Eq. (|H]) corresponding to eigenvalue A in the form 



^2 = ^cCie-®^ + LCse-®!, (21) 

where the parameters Ci and C2 are the arbitrary complex constants, and the other param- 
eters are defined by 

Gi = —-i {kcX — U(t) + D [x + 5t) , 
02 = ~-i {kcX — Uct) — D {x + 6t) , 
L = -^'ikc - D - X, 



D = y(tk, + 2Xf-AAl 
S = -i\ - ^kc, 

Following the formulas (j5]), ( IT8|) . and (|2T|) we can obtain the one- and two-soliton solutions 
embedded in a nonlinear spin wave background, respectively. 

IV. MODULATION OF ONE-SOLITON SOLUTION BY A NONLINEAR SPIN 

WAVE 

Taking the spectral parameter A = Ai = + 11/1/2, here fii and ui are real number, in 
Eq. (1211) and substituting them into Eq. (IT7|) . we obtain the one-soliton solution from Eqs. 
© and (USD as follows 

rrij. = Ac cos (fi + {Q2 simp + Qi cos (fi) , 
rriy = —Ac sin if + {Q2 cos if — Qi sin if) , 



i-(^c + £^gi)'-(|^g2) , (22) 



where 



(f = kcX — Uct, 

61 = 2DiRX + 2 {DiSi)j^t + 2x0, 

$1 = 2Dux + 2{D,6i)jt-2ip^, (23) 
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Li = —ikc/2 — Di — Ai, 
6i = - kj2, 

cos $1, 

Q2 = 2A,Lusinh^i+ (|Li|'-A2)sin<l>i, 

Ai = {\Lif + Al) cosh 01 + 2AcLiR cos ^i, (24) 
xo = -(ln|C2/Ci|)/2, ^o = [arg(C2/Ci)]/2, 

here the subscript R and I represent the real part and image part, respectively. The pa- 
rameters Ci, C2 are the arbitrary complex constants. It should be noted that Ac, kc, xq, (pQ, 
and the complex variable Ai are the adjustable parameters whose values characterize the 
independent solutions. 

When the spin-wave amplitude vanishes, namely A^ = 0, the solution in Eq. fl22]) reduces 
to the solution in the form 

cosht^i 



m, = Jl (25) 

V cosh^^i ^ ^ 

where 

61 = filX + fii^it + 2xo, 

The parameters — 2xo//Ui and 2(/9q/z/i represent the initial center position and initial phase. 
The expression (123]) describes a magnetization precession characterized by the deviations 
/i^ from the ground state, the amplitude a/1 — fif, the velocity — z/i, and the frequency 
{nl -vf) /2-2-2B/ ilSMo) which shows that the magnetic field contribute to precession 
frequency only. In the other hand, when /i^ = 0, the solution fl22|) reduces to a nonlinear 
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spin- wave solution 

rrix = AcCosip, 
iTLy = —Ac sin (yj, 

= (26) 

Thus the solution (122|) describes a soliton solution of ferromagnet embedded in a non- 
linear spin- wave background fl20l) . The properties of envelope soliton is characterized by 
the width 1/ (2Di/j), the wave number of soliton ks = 2Dii, and the envelope velocity 
f 1 = — {Di6i)j^ /Diji. With the expressions of Di and 6i we found that the amplitude Ac 
and wave number kc of spin wave have the different contribution to soliton solutions. The 
width becomes large with the increasing Ac, however, becomes decreasing with the increas- 
ing kc- In the other hand the velocity of envelope soliton increases continuously with the 
increasing Ac, however, the increasing at first, then decreasing rapidly with the increas- 
ing kc- From Eq. fl23|) we can directly see that when Dudu = 5irDir, the parameters 
9i depends only on x which implies the envelope velocity — {Di5i) ^ / Di^i becomes zero, 
i.e., the magnetic soliton is trapped in space by the nonlinear spin wave which shows a 
new technique of management for the soliton. It should be noted this condition is deter- 
mined by the amplitudes /i^, Ac of soliton and spin- wave, and the wave numbers Vi, kc of 
soliton and spin-wave, respectively. It also can be seen that the amplitude of soliton in 
is modulated by the spin wave, and characterized by the spatial and temporal period 
along the line x = — {Di5i) p^t/ Dm — xq/ Dm, denoted by tt (Di^i)^/ [du {Dfj^ + Df^)] and 
ttDir/ [5ii{Djj^ + Djj)] obtained from Eqs. ([22]) and ([MD, respectively. With the 

increasing Ac the envelope of solution hold the shape of soliton, and the envelope valley is 
more deep. Moreover, we can see spin wave move the center position of soliton. 

V. TWO-SOLITON INTERACTION MODULATED BY A NONLINEAR SPIN 

WAVE 

It is well known that soliton has the properties of the physical particles in the process of 
collision. Therefore, to investigate the magnetic soliton collision modulated by a nonlinear 
spin wave is an very interesting phenomenon in spin dynamics. To approach this we should 
obtain the two-soliton solution of m from Eqs. (j7]) and ( ITSj) . Taking the spectral 
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parameter A = = ^j/2-\-iVj/2, j — 1,2, we obtain the two-soliton solution of ferromagnet 
as 

rux = Ac cos </? + ^ {G21 sin (p + G2R cos (/?) , 



ruy = -Ac sin + — {G21 cos (p - G2R sin ip) , 
i'2 



1- Uc + 



\ 2 ^2 



21 



(27) 



where the subscript R and / represent the real part and image part, respectively, and the 
other parameters are defined by 



Oj = 

Ai = 
A2 = 
A3 = 



2D jRX + 2 {DjSj)j^t + 2x0 J, 
2Djix + 2{Dj5j)jt-2^,^^, 

(|Li|^ + A^) cosh 6^1 + 2/1cLircos$i, 
(iLal^ + A^) cosh^a + 2AcL2rCOS $2, 
(L1L2 + ylc) e5(^i~'*i+^2+i*2) _^ g-i(ei-j*i+e2+i*2) 



+^ (Li + L2) 



(28) 



F2 = C1A1A2 - /ii/i2 



A, 



G2 = 

/i = 

/2 = 

/3 = 

/4 = 



/X2C1 A1/2 + /^iCl A2/1 - ^C2/3A3 - ^C2A3/4, 



\2—i^\ 



Ac {L,e'' +L,e-'^) + \L,f e^''^ + A'^e 
Ac {L2e'^ + L2e-'^) + IL2I' e^*^ + Ale-'^\ 

+l2[Lie^(^i+^*i) + A^e-f(''i+i*i)]e-|(«2-i*2)^ 



(29) 
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Ci = |A2 + Ai| , C2 = AtiAt2 (^2 + Ai) , 

D, = sj{ikj2 + \,f -Al 

Lj = —ikc/2 — Dj — Xj, 

5. = -iX.-kj2, (30) 

xo, = -{\n\C2,/C^,\)/2, y^o,, = [arg (C2,,/Ci,,)] /2, 

where Cij, C2J are the arbitrary complex constants, j = 1,2. The subscript R and / repre- 
sent the real part and image part, respectively. It should be noted that kc, Xqj, and 
the complex variable Xj, j = 1,2, are the adjustable parameters whose values characterize 
the independent solutions. In general, the solution fl271) represents the interaction of two 
one-soliton solution modulated by the nonlinear spin- wave (l20i) of ferromagnet. The prop- 
erties of each envelope soliton are characterized by the width 1/ {2Djf>), the wave number 
of soliton 2Djj, and the envelope velocity Vj = — {Dj6j) ^ / Djji, j = 1,2. From Eq. ( !28|) 
we can directly see that when Djj6jj = djuDjn, j = 1,2, the parameters 9j depends only 
on X which implies the velocity of each envelope becomes zero, i.e., the magnetic soliton 
is trapped in space by the nonlinear spin wave. It also can be seen that the condition 
is obviously determined by the amplitudes fij, Ac of soliton and spin- wave, and the wave 
numbers uj, kc of soliton and spin- wave, j = 1,2, respectively. We also observe that the 
amplitude of each envelope soliton in is changed by the spin wave, and characterized 
by the spatial and time period along the line x = — {Dj6j) j^ t/ Djji — Xo/Djji, denoted by 
TT {Dj5j)^ I \bji [p^jR + -D|j)] and 'nDj^j [5^/ {P]r + ■, 3 = 2, for each envelope soli- 
ton, respectively. With the action of spin wave the envelope of solution hold the shape 
of soliton, and the amplitude of each envelope soliton keeps invariability before and after 
collision which shows the elastic collision of two envelope soliton in the nonlinear spin wave 
background. Moreover, we can see spin wave move the center position of soliton oppositely 
before and after collision. 

When /ij = 0, j = 1, 2, the solution fl27j) reduces to the nonlinear wave solution fl20l) . As 
the nonhnear spin-wave amplitude vanishes, Ac = 0, the two soliton solution (127|) reduces 
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to 



{G21 sin7 + G2R cos 7) 
(Gs/ cos 7 - G2R sin 7) 



^2 I 



(31) 



where 
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F2 
G2 



2 1 + 



5 



1 



Ci cosh 6*1 cosh 6*2 — -fJ'ifJ'2 [cos($i — $2) 



cosh(0i + ^2)] 



fi2Ci cosh0i--(C2e'^ + C2e~^^) 



J^2 



+ 



/iiCiCosh^2 - - ((26 



-02 



+ C2e^^) 



The solution fl31l) describes a general scattering process of two solitary waves with differ- 
ent center velocities — z/i and — z/2, different phases $1 — 7 and $2 — 7- Before collision, 
they move towards each other, one with velocity —ui and shape variation frequency Qi = 
inl - ul) /2-2-2B/ (/3Mo) and the other with -z/2 and ^2 = {l4 - ^1) /2-2-25/ (/3Mo). 
Asymptotically, the two-soliton waves flHT]) can be written as a combination of two one-soliton 
solutions (123]) with different amplitudes and phases. The asymptotic form of two-soliton so- 
lution in limits t — )■ — 00 and t — )■ cxd is similar to that of the one-soliton solution (125|) . 
During collision there is no amplitude exchange among three components m^;, my and m^, 
however, a phase and the center position change for each magnetization vector soliton. This 
interaction between two magnetic solitons is called elastic collision. 



VI. CONCLUSION 



In conclusion, we study the dynamics of the magnetic soliton modulated by a nonlinear 
spin- wave in a uniaxial anisotropic ferromagnet. In terms of a reasonable assumption we 
transform the Landau-Lifshitz equation into an equation of the nonlinear type. By means of 
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a straightforward Darboux transformation one- and two-soliton solutions in nonlinear spin- 
wave background are obtained analytically and their properties are discussed in detail. Our 
results show that in the background of a nonlinear spin wave the amplitude of the envelope 
soliton has the spatial and temporal period. The soliton can be trapped only in space which 
is determined by the amplitude and wave number of the magnetic soliton and the nonlinear 
spin wave. The amplitude and wave number of spin wave have the different contribution 
to the width, velocity, and the amplitude of sohton solutions. Moreover, we also observe 
that the envelope of solution hold the shape of soliton, and the amphtude of each envelope 
soliton keeps invariability before and after collision which shows the elastic collision of two 
envelope soliton in a nonlinear spin wave background. 
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